Do Chaotic Trajectories Care About Self-Similarity? 



M. Weiss 1 ' 2 , L. Hufnagcl 1 , and R. Ketzmerick 1 
1 Max-Planck-Institut fur Stromungsforschung and Institut fiir Nichtlineare Dynamik der Universitat Gottingen, 

Bunsenstr. 10, 37073 Gottingen, Germany 
2 EMBL, Meyer hofstr. 1, 69117 Heidelberg, Germany 

We investigate the relation between the chaotic dynamics and the hierarchical phase-space structure 
of generic Hamiltonian systems. We demonstrate that even in ideal situations when the phase space 
is dominated by an exactly self-similar structure, the long-time dynamics is not dominated by this 
structure. This has consequences for the power-law decay of correlations and Poincare recurrences. 
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Generic Hamiltonian systems are neither integrable 
nor chaotic Jjj], but rather exhibit a mixed phase space, 
where regular and chaotic regions coexist. Each island of 
regular motion is surrounded by infinitely many chains 
of smaller islands. As the same holds for any of these 
smaller islands a very complex hierarchical phase-space 
structure is found for generic Hamiltonian systems, which 
is well understood [p| and nowadays appears in textbooks 
on classical mechanics || . 

The dynamical properties, however, are still unclari- 
fied. The most fundamental statistical quantity for char- 
acterizing dynamics is the decay of correlations in time. 
It determines transport properties and is directly related 
to the distribution of Poincare recurrences P(t), which is 
the probability to return to a given region in phase space 
with a recurrence time larger than t. This probability 
decays on average like a power-law M 



P(t) ~ r 



(1) 



due to the trapping of chaotic trajectories in the hierar- 
chically structured vicinity of islands of regular motion. 
The power-law decay is a universal property of Hamil- 
tonian systems. It has dramatic consequences for trans- 
ort H (anomalous diffusion) and quantum mechanics 
(conductance fluctuations and eigenfunctions), which 
sensitively depend on the value of 7. The exponent 7, 
as determined from finite time numerical experiments, 
seems to be non-universal, varies with system and pa- 
rameter, and typically ranges between 1 and 2.5 It 
is a fundamental question of Hamiltonian chaos, how the 
exponent 7 of the dynamics is related to the structure of 
the hierarchical phase space. 

Recently, it was argued by Chirikov and Shepelyan- 
sky that for asymptotically large times the exponent is 
independent of the specific system and parameter and 
is given by the universal exponent 7 = 30. Their 
arguments are based on the universal presence of criti- 
cal tori in phase space and are supported by a numeri- 
cal investigation of the kicked rotor at kicking strength 
K = K c = 0.97163540631. At this parameter value the 
golden torus is critical, i.e., it can be destroyed by an ar- 
bitrarily small perturbation. The self-similar vicinity of 
the critical golden torus [see Fig. |l|(a)] has been studied 



using renormalization methods || and the asymptotic 
value 7 = 3 for the power-law decay of P(t) was pre- 
dicted long time ago ||Jl(J. The fact that it has never 
been observed led to the speculation that the universal 
decay should appear for larger times In Ref. (?J a 

numerical approach allowed to estimate the onset of this 
decay in agreement with the presented data for P(t). 
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Phase-space density of individual trajectories of 
10 7 for the symmetrized standard map. (a,b) 
By stretching the phase space in p direction 
according to the distance p c (q) — p to the critical golden 
torus p c {q) the self-similar phase-space structure in its vicin- 
ity is visualized. The density is determined on a 500 x 500 
grid. (c,d) K — K*: Successive magnifications of the is- 
land-around-island sequence 3 — 8 — 8 — 8 . . .. The density 
is determined on a 250 x 250 grid. In (a) and (c) the trajec- 
tory follows the self-similar phase-space structure, while (b) 
and (d) show representatives of the counter examples. Grey 
shadings are on a logarithmic scale. 



In addition to the sticking of trajectories in the vicin- 
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ity of critical tori, the trapping of trajectories in island- 
around- island structures has been studied (ljjp|. Za- 
slavsky et al. jl3| showed that for the kicked rotor at 
K = K* = 6.908745 the phase space possesses an island- 
around-island structure of sequence 3 — 8 — 8 — 8 . . . [see 
Fig. 0(c)]. They used this self-similarity to derive the 
trapping exponent 7 = 2.25 by renormalization argu- 
ments, which was recently supported by the numerics 

In fact, these renormalization approaches for single 
self-similar phase-space structures can be considered as 
special cases of the more general binary tree model by 
Meiss and Ott fl5|| . In this model a chaotic trajectory 
can at any stage of the tree either go to a boundary cir- 
cle (level scaling) or to the island-around-island struc- 
ture (class scaling) . The universal coexistence of the two 
routes of renormalization at any stage led to the expo- 
nent 7 = 1.96 In contrast, the recent findings claim 
that just one of these scalings is relevant for the trapping 
mechanism: While in Ref. it is argued that universally 
(and in particular for K = K c ) the level scaling should 
dominate, in Ref. jl3| it is claimed that for K = K* the 
class scaling describes the trapping mechanism. 

In order to clarify these contradictions, we numerically 
investigate P(t) for the kicked rotor at K c and K* for 
times larger than in previous studies. We find strong de- 
viations from the predictions of the renormalization the- 
ories that only consider a single self-similar phase-space 
structure. In addition, our numerical approach allows to 
analyze where chaotic trajectories are trapped in phase 
space. For large times the majority of trajectories is not 
trapped in those phase-space regions, that are described 
by the simple renormalization theories. We thereby re- 
veal the mistaken assumption of these theories. In par- 
ticular for K — K c , the self-similar vicinity of the critical 
golden torus does not dominate the trapping mechanism 
for large times and thus even in this ideal situation the 
proposed universal exponent 7 = 3 is not found. For 
K = K* , the majority of long trapped trajectories does 
not follow the self-similar island-around-island structure, 
which leads to a smaller exponent than predicted. Al- 
though the phase space in both cases is dominated by 
exactly self-similar structures, our analysis shows that 
they do not dominate the dynamics. 

We use the standard map (kicked rotor) defined by 

q n +i = q n + Vn mod 2-k p n+1 = p n + K sin q n+1 , (2) 

which has a 27r-periodic phase space in p direction. Wc 
concentrate on two parameters: (i) The dynamics for 
K = K c is bounded in p direction by the golden torus, 
which is critical || . The route towards the critical golden 
torus is determined by the principal resonances given by 
the approximants of the golden mean a = (V5 — l)/2 
and the scaling has been analyzed in detail The dy- 
namics along this route was described by a Markov chain 



leading to 7 = 3.05 Q and alternatively via the scaling 
of the local diffusion rate leading to 7 = 3 jl(| . (ii) The 
phase space for K = K* consists of two small accelera- 
tor modes embedded in an otherwise completely chaotic 
phase space. Each mode shows an island-around-island 
structure of sequence 3 — 8 — 8 — 8 . . . |jl3). This exact scal- 
ing relation was used to predict the exponent 7 = 2.25. 
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FIG. 2. Poincare recurrences P(t) (solid) for the standard 
map with K = K c , deviating from the prediction of Ref. j?J] 
(dashed). The inset shows the decaying fraction f(t) of tra- 
jectories following the route of renormalization towards the 
critical torus. 

In order to check the predictions for P(t), in case (i) 
we start several long trajectories initially located near 
the unstable fixed point (q,p) = (0,0). We measure the 
times r for which an orbit stays close to the critical torus 
by monitoring successive crossings of the line p = 0, as 
was also done in Ref. [Q. In case (ii) we start many trajec- 
tories at four different regions in the chaotic part of phase 
space Whenever a chaotic trajectory is trapped to 
one of the island structures, it follows the dynamics of 
the accelerator mode and jumps to the neighboring unit 
cell in p direction. We measure the time r it continuously 
jumps one unit cell per iteration in the same direction. 
From the set of trapping times r one determines the frac- 
tion P(t) of orbits with r >t. This quantity decays with 
the same power-law exponent as the Poincare recurrences 
P(t) and was chosen for numerical convenience. The total 
computer time corresponds to (i) 15- 10 12 and (ii) 8- 10 12 
iterations of the standard map. We have checked if our 
statistical data for large times are sensitive to the un- 
avoidable finite numerical precision, by comparing data 
for double (wI6 significant digits) and quadruple («32 
digits) precision. We found no difference and present a 
combination of both data sets in Fig. |[ 

In Fig. U we compare our numerical findings for P(t) 
for case (i) with the prediction P(t > 10 7 ) = 3.9- 10 12 i" 3 
extracted from Ref. |?| . This power law is not compatible 
with our data, even though we are in the time regime, 
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where it should be observable according to Ref. TTj . For 
10 5 < t < 10 9 we rather see an exponent 7 « 1.9 Jl7| . 
For case (ii) with K = K* , we find a power-law decay 
of P{t) with 7 = 1.85 (Fig. ||) for various starting con- 
ditions |l6[ c ontradicting the renormalization prediction 
7 = 2.25 |g| . 
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FIG. 3. P (t) for K — K* for various initial condi- 
tions (vertically shifted for better comparison). All 
four curves (solid) show the same power-law behavior in- 
cluding the fluctuations and deviate from the prediction 
7 = 2.25 from Ref. |l^] (dashed). The inset shows the de- 
cay of the fraction f(t) of orbits following the self-similar is- 
land-around-island structure. 

Our numerical results raise the question why the renor- 
malization theories predict the wrong exponents. This is 
particularly surprising since the structure of the phase 
space for the specific parameters K c and K* is dom- 
inated by one exactly self-similar hierarchy. All these 
approaches rely on the assumption that this single self- 
similar phase-space structure also dominates the long- 
time trapping of chaotic trajectories. We will now check 
this assumption. 

To this end we calculate the density in phase space for 
trajectories, that are trapped for long times. For case (i) 
we show two examples in Fig. |l](a,b). Figure |l](a) shows 
a trajectory of length t w 5 • 10 7 that follows the route 
to the critical golden torus up to the principal resonance 
with winding number 55/144. This is consistent with the 
renormalization theory according to the data presented 
in Ref. 0. In contrast, the trajectory shown in Fig. |l](b) 
approaches the critical torus only up to the resonance 
3/8 and is predominantly trapped around a non-principal 
resonance. This trajectory is not captured by the renor- 
malization theory, since it has the same length as the 
trajectory in Fig. ^(a) and therefore should be trapped 
around the 55/144 resonance or one of its neighbors. In 
order to quantify this observation we introduce the frac- 
tion f(t) of trajectories with trapping time t that follow 
the route of renormalization. Numerically, we determine 



f(t) by considering trajectories with trapping times in an 
interval around t. We classify these trajectories, as was 
done in Fig. ||, according to their phase-space densities. 
For t > 2 • 10 this fraction decreases to zero (Fig. ||, in- 
set). This clearly demonstrates that for large times the 
majority of trajectories is not trapped in the self-similar 
phase-space regions approaching the critical golden torus. 
The assumption of the renormalization theories leading 
to 7 ss 3 is violated and thus they are not applicable 
for predicting P(t). From the ratio of the predicted 
P(t) ~ for the trajectories trapped in the self-similar 
phase-space structure and the observed P(t) ~ t~ 19 we 
expect their fraction to decay as fit) ~ t~ 1A for large 
times. This is confirmed by our numerical data in the 
inset of Fig. |[ We find that the majority of trajecto- 
ries is trapped around non-principal resonances, which 
is in agreement with the binary tree model |15| . This is 
in strong contrast to the conclusions of Ref. [J7| that are 
based on the computation of exit times from the vicinity 
of unstable fixed points of principal resonances. As for 
the investigation of local diffusion rates |Q also the anal- 
ysis of the mean exit time yields information only about 
trajectories trapped in the region studied. One cannot 
conclude, however, that these trajectories dominate the 
global trapping mechanism, as is convincingly shown by 
our numerics. 

We carry out the same analysis for case (ii) , K — K* . 
In Figure |l](c,d) we show the phase-space densities of 
two long trajectories. Although both trajectories have 
the same length, only the trajectory in Fig. |l](c) fol- 
lows the self-similar island-around-island structure, while 
the trajectory shown in Fig. |l|(d) is trapped around 
other islands. The fraction fit) of trajectories follow- 
ing the route of renormalization decays to zero for large 
times (Fig. [| inset). This decay is well described by 
the estimate f(t) ~ ^2.25^-1.85 _ f -o.4 [ e ^ the ra _ 

tio of the predicted P(t) ~ t~ 2 - 25 and the observed 
Pit) ~ i -185 . This shows that the renormalization the- 
ory for the island-around-island structure is not capable 
of explaining P(t). It should be noted that this difference 
is not caused by the effect, that the finite precision of K* 
eventually leads to a breakdown of the self-similarity on 
very small scales. 

In conclusion, our analysis shows that even in the pres- 
ence of an exactly self-similar phase-space structure it 
is not sufficient to describe the trapping mechanism of 
chaotic trajectories by only this structure, as was re- 
cently claimed in the literature. We find that additional 
island structures may dominate the trapping mechanism 
for large times and thus affect the power-law decay of 
P(t). Our analysis supports qualitatively the tree model 
by Meiss and Ott, which allows for the coexistence of two 
routes of renormalization at any stage. Quantitatively, 
we find slightly smaller exponents. It remains an open 
question if there exists a universal asymptotic exponent 
for the trapping of chaotic trajectories in Hamiltonian 
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systems. 
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